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Abstract 

As a consequence of the large top quark Yukawa coupling, supersymmetric unified theories 
with soft supersymmetry breaking terms generated at the Planck scale predict lepton flavour 
and CP violating processes with significant rates. 

The flavour violating parameters of the low energy theory are derived in both SU(5) 
and SO(IO) theories, and are used to calculate the rate for ^ —>■ ej. The sensitivity of 
the search for /i — > 67 is compared with that for /j, —>■ e conversion in atoms, r — > /i7 and 
the electric dipole moment of the electron. The experimental search for these processes is 
shown to provide a very significant test of supersymmetric unification, especially in SO(IO) 
but also in SU(5). 
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1 Introduction 



The importance of looking for direct tests of unified theories cannot be overestated. As is well known, 
such an opportunity is essentially restricted to the study of violations of those conservation laws which are 
valid in the Standard Model as a consequence of exact "accidental" global symmetries. We refer to baryon 
number, B, and to the individual lepton numbers, Le, and Lr- 

In these respects, the violation of individual lepton numbers while preserving the overall lepton number, 
L = Le + + Lt, — hereafter called Lepton Flavour Violation (LFV) — plays a special role. If the Grand 
Unified Theory, characterized by a large mass scale Mq, has the pure Standard Model as its low energy 
approximation, the rates for the corresponding LFV processes (/x ej, /j, ^ e conversion, ^ 3e, etc.) 
are unobservably small, since they are necessarily mediated by non-renormalizable effective interactions 
scaled by inverse powers of Mq. On the contrary, in a supersymmetric unified theory with supersymmetry 

— 1/2 

effectively broken at the Fermi scale, m = 0{Gp ), the rates for the LFV low energy processes are 
only suppressed by powers of 1/m Q. In general this would actually also be the case for B and/or 
i-violating processes, like proton decay, strongly suggesting the need of matter parity (or it!-parity) in a 
unified supersymmetric theory. Correspondingly, the LFV processes, consistent with matter parity unlike B 
and/or L violations, emerge as very interesting possible experimental signals of supersymmetric unification. 

In a previous paper Ql, two of us (R.B. and L.H.) have pointed out that the large Yukawa coupling of 
the top quark at the unification scale, Xtc is an important source of flavour violation which reflects itself, 
via the unified couplings, in relatively large rates for general LFV processes. There it has been argued 
that the study of the corresponding experimental signals provides a test of supersymmetric unification at 
least as significant as the one that can be obtained from either proton decay or neutrino masses. In the 
present work we substantiate further this statement, by making an analytic study of the rates for one of the 
processes discussed in ref. [||, ^ — > 67, in the full parameter space of the unified theory, both in SU(5) 1^ 
and in SO(IO) [Q and by subsequently comparing it with the other processes and quantities of interest. In 
addition to generating lepton flavour violating interactions, the large top quark Yukawa coupling also leads 
to important contributions to the electric dipole moment of the electron and neutron in SO(IO) theories. 
In this paper we give the electric dipole moment over the full parameter space of the SO(IO) unified theory. 

The paper is organized as follows. In section 2 we discuss the physical mechanism which allows the top 
quark Yukawa coupling to generate large amplitudes for processes which violate individual lepton numbers. 
In section 3 we summarize the present information on Ajg and we describe an upper bound on A^g arising 
from the infrared fixed point behaviour of the top Yukawa coupling above the unification scale up to Mpi. 
In section 4 we study the scaling of the supersymmetry breaking parameters in SU(5) with emphasis on the 
flavour violating effects due to Xtc- In section 5 we give, in SU(5), the pieces of the low energy Lagrangian 
relevant to the calculation of the LFV processes in the physical lepton and slepton basis. In section 6 we 
calculate the rate for /i — > 67 in the full space of parameters. In sections 7-9 we extend the analysis of 
sections 4-6 to the SO (10) case. In section 10 we discuss the fi ^ e conversion in atoms and we establish 
the relative merit of the study of this process with respect to /x ^ 67 in the search for a signal of lepton 
flavour violation. The same is done in section 11 for the r — + fij decay. Finally, in section 12, we study the 
relation between n ^ and the electric dipole moment of the electron |^ . Our conclusions are drawn in 
section 13. Appendices A and B contain the analytic solutions of all the relevant Renormalization Group 
Equations from Mp\ to Mq (appendix A) and from Mq to Mz (appendix B), both in SU(5) and in SO(IO). 

2 The origin of lepton flavour violation 

In this paper we study grand unified theories which incorporate weak-scale supersymmetry j6| and have the 
origin of supersymmetry breaking near at the Planck scale 0. These theories lead to the successful weak 
mixing angle prediction, and, as a promising direction for unifying both the forces and the fundamental 
fermions, are currently receiving much attention. In all such theories, we find that the large top quark 
Yukawa coupling leads to a rate for fj, which can be reliably computed in terms of weak-scale 

parameters |2] . Over much of the interesting parameter space, the rate is within two orders of magnitude 
of the present experimental limit. At first sight, it is surprising that the top quark Yukawa coupling should 
lead to any violation of Le or L^. What is the physical origin of this effect, and why is it not suppressed 
by inverse powers of Mq7 The answer lies in new flavour mixing matrices, which are analogous to the 
Cabibbo-Kobayashi-Maskawa matrix. 



In the standard model the quark mass eigenstate basis is reached by making independent rotations on 
the left-handed up and down type quarks, and d^. However, these states are unified into a doublet of the 
weak SU(2) gauge group: Q = {u^.d^). A relative rotation between ul and therefore leads to flavour 
mixing at the charged W gauge vertex. This is the well-known Cabibbo-Kobayashi-Maskawa mixing. With 
massless neutrinos, the standard model has no analogous flavour mixing amongst the leptons: the charged 
lepton mass eigenstate basis can be reached by a rotation of the entire lepton doublet L — {v^, e^). 

How are these considerations of flavour mixing altered in supersymmetric unified theories? There 
are two new crucial ingredients. The first is provided by weak-scale supersymmetry, which implies that 
the quarks and leptons have scalar partners. The mass eigenstate basis for these squarks and sleptons 
requires additional flavour rotations. As an example, consider softly broken supersymmetric QED with 
three generations of charged leptons. There are three arbitrary mass matrices, one for the charged leptons, 
e, and one each for the left-handed and right-handed sleptons, cl and eu. To reach the mass basis therefore 
requires relative rotations between cl and ei, as well as between cr and cr, resulting in two flavour mixing 
matrices at the photino gauge vertex. 

In supersymmetric extensions of the standard model, these additional flavour-changing effects are known 
to be problematic. With a mixing angle comparable to the Cabibbo angle, a branching ratio for — > 67 of 
order 10"'* results. In the majority of supersymmetric models which have been constructed, such flavour- 
changing effects have been suppressed by assuming that the origin of supersymmetry breaking is flavour 
blind. In this case the slepton mass matrix is proportional to the unit matrix. The lepton mass matrix 
can then be diagonalized by identical rotations on and ei, as well as on cr and cr, without introducing 
flavour violating mixing matrices at the gaugino vertices. Slepton degeneracy renders lepton flavour mixing 
matrices non-physical. 

The unification of quarks and leptons into larger multiplets provides the second crucial new feature in 
the origin of flavour mixing The weak unification of iti, and dL into Q is extended in SU{5) to the 
unification of Q with and e£ into a 10 dimensional multiplet T^Q^ul, el). Since higher unification 
leads to fewer multiplets, there are fewer rotations which can be made without generating flavour mixing 
matrices. 

In any supersymmetric unified model there must be at least two coupling matrices, Ai and A2, which 
describe quark masses. If there is only one such matrix, it can always be diagonalized without introducing 
quark mixing. One of these coupling matrices, which we take to be Ai, must contain the large coupling, 
At, which is responsible for the top quark mass. We choose to work in a basis in which Ai is diagonal. 
The particles which interact via Af are those which lie in the same unified multiplet with the top. In all 
unified models this includes a right-handed charged lepton, which we call . This cannot be identified as 
the mass eigenstate r£ , because significant contributions to the charged lepton masses must come from the 
matrix A2, which is not diagonal. 

The assumption that the supersymmetry breaking mechanism is flavour blind leads to mass matrices 
for both cl and cr which are proportional to the unit matrix at the Planck scale, Mpi. As we have seen, 
without unified interactions, lepton superfield rotations can diagonalize the lepton mass matrix without 
introducing flavour mixing matrices. However, the unification prevents such rotations: the leptons are in 
the same multiplets as quarks, and the basis has already been chosen to diagonalize Ai. As the theory 
is renormalization group scaled to lower energies, the At interaction induces radiative corrections which 
suppress the mass of cr^ beneath that of br^ and br-^ . Beneath Mq the superheavy particles of the theory 
can be decoupled, leaving only the interactions of the minimal supersymmetric standard model. Now that 
the unified symmetry which relates quarks to leptons is broken, a lepton mass basis can be chosen by 
rotating lepton fields relative to quark fields. However, at these lower energies the sleptons are no longer 
degenerate, so that these rotations do induce lepton fiavour mixing angles. Radiative corrections induced 
by At lead to slepton non- degeneracies, which render the lepton mixing angles physical |^]. 

This discussion provides the essence of the physics mechanism for lepton fiavour violation in superunified 
models. Since the flavour mixing matrices have complex entries, they also lead to CP violation. It shows 
the effect to be generic to the idea of quark-lepton unification, requiring only that the superpartners have 
masses around the Fermi scale, and that supersymmetry breaking be present at the Planck scale. The 
imprint of the unified interactions is made on the soft supersymmetry breaking coefficients, including the 
scalar trilinears, which are taken to be flavour blind at the Planck scale. Eventually this imprint will be 
seen directly by studying the superpartner spectrum, but it can also be probed now by searching for L^, 
L^, Lr and CP violating effects. 



3 The top Yukawa coupling at the GUT scale 



The top Yukawa coupling at the unification scale, AtCi plays a crucial role in the determination of the LFV 
effects discussed in this paper. In this section we therefore summarize the present information on Ajg which 
comes from two different sources: the direct measurement of the top mass and, indirectly, the bottom/tau 
mass ratio. 

The top Yukawa coupling Xtc can of course be easily scaled down to determine its value at the weak 
scale At (see eq. (54b) of app. In turn, At determines the top quark pole mass via ||] 



where, as usual, tan/3 = Vu/v^ is the ratio of the two light Higgs vacuum expectation values. Figure |i| 
shows AtG as function of the strong coupling constant az{Mz), for Mt = 174 ± 16 GeV and for moderate 
(tan/3 = 2) or relatively high values (tan/3 = 10) of the ratio Vu/wd- The rapid saturation for large tan/3 
implies that the lowest curve in fig. |^ is actually a lower bound on Ate for mt > 158 GeV. 

As is well known, the Yukawa supcrpotential of minimal SU(5) allows a prediction for the ratio mh/m^ 
as a function of Afc and a^{Mz) 1^, This prediction is given for the running 6-mass nrn,{mi,) in fig. |^, 
and is compared with the preferred value as determined from T-physics. The dependence on tan [3 drops 
out in the ratio mj^/m^^ except for tan/3 w mt/mi,. For moderate values of tan /3 there is a clear consistency 
between fig. ^ and fig. |^ with a strong indication for a high value of Ate- The consistency is weaker for 
larger tan /3, unless the top mass is close to 200 GeV, in the upper range of preliminary values indicated by 
the CDF experiment ||ll|. In this case, of course, a rather high value of AtG is also indicated, resulting in 
a large fiavour violation of the soft supersymmetry breaking parameters at the unification scale. For very 
large tan /3, close to mt/mb, consistency with mb/nir is possible even for the smallest values of Ate allowed 
by the top mass. We have not studied this case in this paper. However, because the fi ej amplitude 
always contains a term proportional to tan/3, the rate is always significant for such large values of tan/3. 

For later purposes, it will be useful to know the behaviour of At at energies above Mq, as determined 
from the Renormalization Group Equations (RGE). Assuming that the unified gauge coupling and At 
itself are the only relevant couplings, the RGEs are solved in appendix A at the one loop level. The solution 
for At displays an infrared fixed point. If At(A/pi) is large, but still perturbative, it will be drawn to the 
infrared fixed point at Mq. The value of the coupling at the fixed point at Mq is larger for SU(5) than 
for SO(IO), and depends on the one loop coefficient, Bq, of the gauge /3-function, as shown in figure I The 
quantity Xf^^ plotted in this figure is the value of Ate for which the one loop evolved value of At(Afpi) 
becomes infinite. For all the numerical work of this paper we take Ate < •^tcf'^j so that perturbation theory 
can be trusted. For larger values of Ate the theory becomes non-perturbative at scales beneath Mpi. 
Although we are unable to make computations for this case, the non-perturbative coupling is expected to 
generate large non-degeneracies amongst the scalars, leading to large rates for /i — > ej. 

4 Scaling of supersymmetry breaking parameters in SU(5) 

The messengers of flavour violation in the Icpton sector are the soft supersymmetry breaking terms, which 
are therefore crucial to determine. Without having to specify the actual mechanism of supersymmetry 
breaking, nor the sector in which it takes place, we assume that it is transmitted to standard matter by 
supergravity couplings and that it results, at the Planck scale, in universal soft breaking terms. 

Standard matter occurs in the usual triplication of 10 (T) ® 5 (F) representations of SU(5), which are 
coupled to a 5 (H) and a 5 (H) representation of Higgs supermultiplets in the Yukawa supcrpotential 



which we assume to be valid from Mq to the Planck scale. The full supcrpotential will contain other 
supermultiplets S, needed to break SU(5) but not directly coupled to matter. Assuming no large Yukawa 
couplings of the S fields to the H, H multiplets, the S fields affect the determination of the soft supersym- 
metry breaking terms at the GUT scale only via their contribution to the gauge /3-function from Mq to Mp\. 
Unless otherwise specified, we shall take the SU(5) /3-function coefficient of the minimal Dimopoulos-Georgi 
model H]. Different /3-function coefficients mostly affect the rates for the LFV processes only through the 
restrictions that they induce on the range of the low energy parameters (see appendix ^ . 




V ^ 174 GeV, 



(1) 



(2) 



Figure 1: The top Yukawa coupling at Mq for tan/3 = 2 (full lines) and tan/3 — 10 (dashed lines) for 
Mt = 158, 174, 192 GeV (in increasing order), as function of a^{Mz)- 



Figure 2: The running 6-quark mass in the a^{Mz),\tG plane from h/r unification. The darker area 
corresponds to ■mt,{mi,) ~ 4.25 ± 0.10 GeV, as obtained from T-physics. 

Figure 3: Fixed point upper bounds on Ate in SU(5) and SO(IO), as defined in the text, as functions of 
the one loop coefficient bQ of the gauge /3-function. 

The relevant part of the soft supersymmetry breaking Lagrangian, before SU(5) breaking, has the form 

- Csott = Vsoit = f tm|f +F^m?pF + m]j\H\'' + m\\H\'' + f^A'^X'TH + f^A^X^FH (3) 
with, at the Planck scale, 

= m?p = itiqI, 

= "ifi = "^o> (4) 
A" = A'* = Aol. 

The renormalization of the parameters in (|^) down to the GUT scale is most easily done by working in 
the basis where the Yukawa matrix A" has diagonal form (hereafter called the u-basis). By keeping in the 
RGE only the one loop effects due to the SU(5) gauge coupling and to the third entry of A", A33 = At, it 
is simple to rescale down to Mq the soft breaking parameters (see appendix ^) . Flavour universality is of 
course no longer maintained. In fact, the mass term for the ten-plet and the A-terms acquire the form 

m^Q = diag(m|G, m^Q, m|(3 - Iq) = vo^qI - /q, (5a) 

A% = diag(ylrfG, ^dG, ^dG - \Ig) = ^^gI - ^Iq^ (5b) 

= diag(A„G - ^Ig^^uG - j^g- - Iq)- (5c) 
At the same time, the mass matrix of the five-plets maintains the universal form 

= mpcl (6) 
and the Yukawa coupling matrix A^ remains diagonal, 

Aq = diag(A„G , Acg , Af g ) , (7) 

whereas A'* gets renormalized to Aq- The explicit expressions for tti^q, rripQ, AdG, ^uG, Ig and Iq, as 
well as the renormalization of the Higgs mass parameters are given in appendix The flavour breaking 
parameters Iq and I'q have a crucial dependence on the top Yukawa coupling at Mg, XtG (see section 5). 

In this paper we take a universal boundary condition for the soft supersymmetry breaking parameters 
at the Planck scale. How do our results depend on this assumption? A well motivated relaxation of this 
assumption is to allow soft scalar masses to be the most general allowed by the gauge symmetry and by 
a symmetry which interchanges one generation with another. This would satisfy flavour changing phe- 
nomenology without forcing identical Higgs and matter scalar masses, and would also allow the scalars in 
T to have masses different from those in F. Although extra parameters must be introduced, this general- 
ization will not affect our results in a crucial way. More important would be the addition of small flavour 
changing scalar masses at the Planck scale, since they would lead directly to the processes which we discuss 
in this paper. These contributions would simply add to those which we calculate here. While cancellations 
cannot be excluded, we believe they would have to be accidental. For example, the contributions from the 
Planck scale boundary condition would arise from string physics and would be independent of the value of 
Mg- On the other hand, the contributions calculated in this paper do depend on Mq. 

5 The low energy Lagrangian in SU(5) 

After SU(5) breaking, the scaling down to low energy of the various parameters results in the low energy 
Lagrangian, whose relevant pieces are summarized for ease of the reader. They are, to first order in the 
Yukawa couplings: 



i. The slepton mass matrix 



- £^ = Vm\L + e^j^mleR + e^iA" + l^itan [3)X'eLVd + h.c. (8) 

where L, cr are 3- vectors containing the SU(2) doublet and singlet sleptons, 

ml ^ mil, ml^mll-lG, A^^A.l-^I'^, (9) 

ml, m\ and A^. are given in appendix and a term proportional to the /i parameter has been 
explicitly introduced; 

ii. the Higgs mass terms 

-Cm = {ml + + [ml + Ai^)|/idP - ml^{Khd + h.c.) (10) 

with m^, given in eq. (62) of appendix^. 

iii. the quarks and lepton mass terms 

Cy - Q^\lul ■ Vu + Q^\%dl ■ Vd + efX'zL ■ (11) 

where, in the u-basis, has kept its diagonal form and the matrices A'^ and A°, equal at Mq, have 
been shifted by the different renormalization effects due to At and the gauge couplings. 

The LFV parameter Iq is directly related to the splitting between the tr and the cr (/ifl). The ffl-mass 
is shown in figures ^ for fixed values of Ajg and of the e^-mass, as function of and of the wino mass M2 
in its full range, as determined from mg^ itself. We take the value of Ajg such that A^q = 0.8{X^^)^ . The 
lightness of tr is of course a main consequence of the present picture as far as the superpartner spectrum 
is concerned. Another interesting consequence is the strong upper bound on the gaugino mass for any 
given mg^, which results in particular in the lightest supersymmetric particle being always the lightest 
neutralino. The -r^-mass has a negligible dependence on tan/3 for to|^ ^ . For values of mg^^ higher 
than 300 GeV, m^^ and M2 rescale in the same way as mg^ itself does. 
By diagonalizing A^ and A^, we have 

A|wd = V*M'^U^ (12a) 
A|t)d = F°*M''[7'=t (12b) 

where Af^, M° are the diagonal mass matrices for down quarks and charged leptons, U = U°, V is 
the usual Cabibbo-Kobayashi-Maskawa matrix and, as an effect of the top Yukawa coupling, the matrix 
elements of V° are related to those of V by |l^ 

V^j — yVij for i ^ j and {i or j) — 3, V^j — Vij otherwise (13) 

and y is defined in eq. (|5^ ) of app. B. We ignore for the time being the fact that one does not obtain in this 
way the correct relation between the masses of the light leptons and down quarks, which are also related 
to each other by an appropriate renormalization group rescaling. 

It is convenient to work in a mass eigenstate basis for the charge leptons, which is simply obtained by 
the redefinitions (with primed indices suppressed after eq. ([l^)) 

In the gaugino couplings, the rotation on the charged lepton doublets can be compensated by the same 
rotation of the full supermultiplets, since the SU(2) doublet slepton mass matrix has kept its diagonal form 
through renormalization, whereas this is not the case for the singlets er. As a consequence, the matrix V 
appears in the neutralino couplings 

Cg = ^g'Y, [- -erei7V„(iJ„^ +cot%i7„,^J + 4y^tg^^^^^^ (15) 



Figure 4: Isoplots of mf^, in the M2, Ae/rrii^ plane for (a) Ate — 1-4, mg^ — 100 GeV, tan/3 = 2 in SU(5), 
(b) AtG = 1-4, mg^ = 100 GeV, tan/3 = 10 in SU(5), (c) Ate = 1-4, mg„ = 300 GeV, tan/3 = 2,10 in 
SU(5), or, for M2 < 270 GeV, Ate = 1-25, mg„ = 300 GeV, tan/3 = 2,10 in SO(IO), (d) Ate = 0.86, 
niej^ = 300 GeV, tan/3 = 2, 10 in SO(IO). In fig.s ^,b the isohnes are separated by 10 GeV, in fig.s ||c,d 
by 30 GeV. 

Figure 5: Lepton flavour violating couplings in SU(5). 

Figure 6: Diagrams giving rise to the decay fi —> ej in SU(5). In figures ^,c, an external photon line is 
left understood, which can be attached to either of the scalar lines. 

Figure 7: Isoplots of B.R.(/i ej) in SU(5) in the M2, Ae/me„ plane for Xto ~ 1-4, mg^ — 100 GeV and 
(a) tan/3 = 2, fi < 0, (b) tan/3 = 2, ^ > 0, (c) tan/3 = 10, ^ < 0, (d) tan/3 = 10, ^ > 0. The dashed 
(dotted) lines delimit regions where m?^ < (/i^ < 0). The shaded area also extends to m^^ < 45 GeV. 
The darker area shows a region where the rate is small, and passes through zero, due to a cancellation of 
terms. The dot-dashed line corresponds to the present experimental limit. For the CKM matrix elements 
we take \Vab\ = 0.04 and \Vtd\ = 0.01. 

Figure 8: Same as in fig. ^ for mg^ — 300 GeV. 



where Nn are the four neutralino mass eigenstates, of mass M„, related to the bino and the neutral wino 

by 

Notice that, in the slepton basis in which we are working, also the third term in the right side of (^) has 
non diagonal form, being 

- C'f^'^- = {Ae + //tan/3)g]Jy"M'=gL - 4i7^V"*M^gL + h.c. (17) 



6 /i ^ 67 in supersymmetric SU(5) 

The LFV couplings are summarized in fig. ^. Correspondingly, if we neglect the electron mass and we 
work to first order in TOg/rn^ or nifj^/mr, the diagrams giving rise to the decay — > 67 are shown in fig. |^. 
Taking into account that the selectron, e^j, and smuon, flu, singlets are degenerate at a common squared 
mass m| , whereas they are split from the stau singlet f^, of squared mass m?^ = rn|^ — Iq, and using 
the unitarity of the matrix V, one obtains the following contributions to the fi ej decay amplitude 

f — T5 

A^^ifi ^ e^) = -ie- Ueicr^^q'' — - — (18) 



,(a) 



2 - 4^cos2 0w 
b. from the diagram of figure 



m^V:;v:,[G,{mlJ ~ G^imjJ] (19a) 



4:71 cos^ 6 



w 



m,Vr^V,^Me + ^^tanp)[G2{mj^,mlJ - G^imj^^mjj] (19b) 



c. from the diagram of figure ||c; 

a 

An cos^ 6w 



^2"^ = .,.^.2.,.. "^MK"K''e(-i^^)G2«,m?J (19c) 



where 



and 



n=l ^^^^ 

r. I 2 2n G2{ml)~G2{ml) I 2,011 



'1 '"2 

Correspondingly, the decay rate is given byll 



2(r- l)s 



r(M - 67) - IF2 r , F2 = F!t^ + ) + . (20) 

Equations (19j20|), together with the expressions of the parameters in the low energy Lagrangian as defined 
in the previous section and explicitly given in appendix allow the numerical calculation of the branching 
ratio B.R.(/i — > 67), shown in fig. ^and fig. ||for mg^ equal to 100 GeV and 300 GeV respectively. For values 



of rnep. greater than 300 GeV and fixed Ag/mg^, Milm^j^ the branching ratio scale as m~ . From ( |l3[) , 



for the CKM matrix elements we take |Vcfc| = 0.04 and |Vtd| = 0.01. 

The set of independent parameters, on which the branching ratio depends, are {Ao,mQ,M5Pi} which 
determine the soft operators, the top quark Yukawa coupling, the coefficient of the one-loop gauge beta 
function of the unified theory, 6g, the ratio of weak vacuum expectation values, tan /3, and the Higgs mixing 
parameter \x. We choose to exchange {Aq, ttiq, M5P1} for the physically more interesting set {Ag, to|^, M2}, 
where is the light generation lepton A-parameter, mg^j is the mass of the right-handed selectron and M2 
is the weak scale gaugino mass parameter for SU(2). In appendix A the full dependence of quantities on 
fee is given, and is found to be mild, hence we have chosen the minimal value 6g = ~3. The /i parameter 
which enters in ^"2*^^ ^Q- (19t)), and also in the neutralino mass matrix, is expressed, up to its sign, in terms 
of the other parameters by means of the electroweak symmetry breaking relation 

2__A/| ra\ - ml tan^ /3 
^ " 2 l-tan2/? ^^'^ 

with m^, TO^ given in appendix H. 



7 Scaling of supersymmetry breaking parameters in SO (10) 

In the case of SO(IO) gauge symmetry, in fact as one of its most attractive features, the quarks and leptons 
of a single generation are the components of a single 16-dimensional spinorial representation This is a 
crucial feature for the problem at hand; it causes all the scalars of the third generation, and not only those 
in the 10 of SU(5), to be lighter than the corresponding scalars in the first and the second generation. In 
turn, and at variance with SU(5), LFV interactions arise also involving the left handed sleptons. With this 
in mind, the considerations of the previous sections can be straightforwardly extended to the SO(IO) case, 
after specifying the Yukawa superpotential and the gauge /^-function, at one loop, at the Planck scale. For 
simplicity we assume that SO(IO) is broken at once to the low energy standard group at Mq. 

In SO(IO) gauge theories a single Yukawa interaction of the three spinorial matter multiplets '^i to a 
vector lO-dimensional Higgs representation $, ^'-^A5'<I>, does not describe any intergenerational mixing, 
since 4* can be rotated to make A diagonal. To describe the mixing, we introduce two 10-plets, and $d, 
in the superpotential Q 

W'so(io) = *^A"*$u + I'^A'^l'^d (22) 

and we assume that the light Higgs doublets (with weak hypercharge Y — -t-i/z) and /id (with weak 
hypercharge Y — — V2) lie respectively in <I>u and $d- As in the SU(5) case, this superpotential is taken to 
be valid already at the Planck scale. Furthermore, here too it is preferable to work in the basis where A", 
which is responsible of the Q — 2/3 quark masses, is diagonal. In analogy with equations (^, at the Planck 
scale we take 

m% = mil, = ml^ = mg. A" = A'' = AqI. 

^In the limit of small Iq, this rate agrees with the analytic expression given in rcf. However, in view of the values 
actually taken by 7q, the expansion generally gives a poor approximation to the correct rate. Previous calculations of the 
^ 67 rate for special values of the gaugino and slepton masses were made in references [h3| . 



After renormalization at the unification scale, we have 



- Aoft = ■t^soft = ^WIg^- + m|^G|$,|2 + ml^cl^dP + ^^A^A^^-^u + ^^A^A^^^d (23) 

where 

171%^ = diag(m|Q,m|Q,m|Q - /g) = mlgl - /g, (24a) 

A^ = diag(ArfG,AdG,^dG-f/G) = ^<iGl-f/G, (24b) 
A^i; = diag(A„G"f/G'^-G-f/G'^-G-/G). (24c) 

When possible we also keep the same notation as in the SU(5) case, but of course the relations of the 
various quantities, e.g. Iq in eq.s (24a) and (5a), to the input parameters at the Planck scale change. 
These relations in the SO(IO) case are given in appendix as function of the one loop coefficient of the 
gauge /3-function. In the text we take bo = —3. 



8 The low energy lagrangian in SO (10) 

The further scaling down from Mq to the weak scale of the different parameters gives rise to the low energy 
lagrangian with the same form as in eq. (|^) , except that now also the diagonal squared mass matrix of the 
left handed sleptons has a split third eigenvalue 



'■LG 



(25) 



The masses of the third generation sfermions are all reduced relative to the ones of the first two generations. 
For example the ffl-mass shows approximately the same pattern as in the SU(5) case for a correspondingly 
lower value of AtG by a relative amount A5^^'^(SO(10))/A5^^'^(SU(5)) « 0.87 (see fig. |). 

On the other hand, in the ferniion mass terms, the symmetry in flavour space of the SO(IO) coupling 
16i 16j 10 gives rise to a symmetric lepton (or down) mass matrix, so that, in eq. (12), If^ = V. 

As before, to calculate the amplitudes for the LFV processes, it is convenient to go to the mass-eigenstate 
basis for the charged leptons. At variance with the SU(5) case, however, this time the term Iq in ( p5| ) 
prevents a counter-rotation also in the left-handed sleptons. As a consequence the flavour changing matrix 
V appears in all the gaugino couplings, which acquire the form (for all terms involving the charged leptons) 



4 



- ^eZV^CLNniH^g + cot OwH^^wJ + 4 V'^inNnH^B + ^i-C- 
ei:V°^i>LiXcK^^)+h.c. 



where Xc are the two chargino mass eigenstates, related to the charged wino by 



cW 



(26) 



(27) 



and i>L is the 3- vector of the left-handed sneutrinos, which, apart from SU(2)(g)U(l) breaking, are degenerate 
with the charged left-handed sleptons. 

Finally, as in SU(5), there is still a non diagonal 'chirality breaking' scalar mass term 



- r^'^- = (^e + A* tan/3)e^y"M"y°tg^ _ i^lf/^, M'^V^jiL + h.c. 
All the LFV couplings in the SO(IO) case are summarized in fig. ^. 



(28) 



9 /i ^ 67 in supersymmetric SO (10) 

The Feynman diagrams contributing to /i — s- ej in SO (10) are shown in figure |l^ (for vanishing electron 
mass). The ones in figures |lO|afl,a^,d, with the helicity flip in the external fermion line, give an amplitude 
proportional to the muon mass, whereas the diagrams of flgures |l^L,ij,Ci.fl,c^ ^, with the helicity flip on 
the internal fermion line, have a dominant term proportional to the tau mass. As such, they dominate the 
decay rate in all of the physically allowed space of parameters. In the approximation of only keeping the 



Figure 9: Lepton flavour violating couplings in SO (10). 

Figure 10: Diagrams giving rise to the decay /i — > 67 in SO(IO). The graphs [lO|a^,b^,CL,c^ involving an 
external right handed muon and an internal neutralino are not displayed. As in fig. |^,c, in fig. |l^,c,c' an 
external photon line is understood. 

Figure 11: Isoplots of B.R.(/i ej) in SO(IO) for me„ = 300 GeV, Ate = 1-25 and all other parameters as 

in fig. g. 

Figure 12: Same as in fig. |ll| for AtG = 0.86. 
Figure 13: Isoplots of C.R.(/i ~> e in Ti) in SU(5) for mg„ = 100 or 300 GeV, Xtc = 1-4 and tan/3 2. 



Figure 14: Box diagrams contributing to 11 —f e conversion in SU(5). 



terms proportional to rrir, the left-handed and the right-handed muon have equal decay amplitudes, from 
the diagrams [lo|bfl,Cfl,Cp and P^L;CL,Cr respectively, which however do not interfere with each other for 
vanishing electron mass. 

From the diagram of figure one has 



^2'"^ = — -^m.K;K=e(K")'(^e +/^tan/3)x 

47rcos"^t^w (29) 

X [G2 (m?^ ,mjj- G2 {mj^ ,mjj~ G2 (m?^ ,mjj + G2 {mj^ , mjj] , 
whereas, from the diagram of figure |l0|cji',Cp one has 

4tt cos^ t/w (^30) 
X [G2(mi„m?J - ^G2{ml,mlJ - iG2(™|„miJ]. 

For the decay rate one has 



a 



r(M-e7) = ^m3|i^2p, ^ F^'"'^ + F^''^^ + Ff'^l (31) 

The isoplots of B.R.(^ 67) are shown in figures |ll|, |l^. 

10 II ^ e conversion 

Keeping only the vector coupling to the nucleus A/", the general amplitude for fi ^ e conversion process 
can be written as 

A = ie^iATj^^ [ueigiRl^^VR - 92r'-^^Vl)u^] + (i? ^ L). (32) 
This amplitude gives rise to the coherent conversion rate 

Tifi ^ e) = Aa'^\F{q)\'ml{\g^R - 52flP + \9il - 92l\'). (33) 

where Z is the charge of the nucleus, Zcs is an efi'ective charge and F(q) the nuclear form factor p^ . 

In our case, the amplitude receives contributions both from penguin-type (P) and from box (B) diagrams. 
More precisely, to leading order in the lepton and light quark masses, the penguin diagrams contribute both 
to gi and 92, unlike the box diagrams, which only contribute to 91 

51 =.9? +5^, 52=52- (34) 
If we define, in analogy with eq. (|l8|), the general off-shell ^ —> amplitude .4p(/i 67) as 

Af_,{n ~* 67) = -ie ■ Ue[q^ FiR-f^J^R + ia^j.^q" F2rVl\u^_, + {R ^ L) (35) 



one has 



9lL,R = ZFiL,R. 



p 

92L,R 



ZF2L,R/'m^. 



(36a) 
(36b) 



In the case of the SO (10) gauge theory, only the magnetic penguin-type amphtudes g^^ ^ have a term 
proportional to nir, as discussed in section 9. Furthermore g^^ ~ g^j^ = 32 • Therefore, apart from terms 
of relative order (m^/mr)^, in view of eq. (36b), 



Tifi ^ e) = 16a^Z^f,Z\Fiq)\'T{t, ^ ej). 



(37) 



In the case of Ti|f, for which jT4| Z = 22, Zcs — 17.6, |-F(9)| = 0.54, taking into account the experimental 
value for the capture rate r{fi capture in Ti) = (2.590 ± 0.012) • 10^/sec we obtain 

C.R.if^ ^ e in Ti) Ti) ^ ^ ^ ^^ ^ ^^^^ 

1 (/i capture m li) 

or, normalizing both ratios to the present upper limits Q 

C.R.(/x ^ e in Ti) BJlV^fT) 

iFT^ = ^-^^ 4.9. 10-" ■ 

Using relation (^ ) , the contours of figures |l]and |l| can be relabelled with values of C.R.(/i — > e in Ti). 
This relation holds whenever the gi form factor contribution can be neglected. In the SO(IO) model this 
is always the case, giving this relation a wide applicability; wider than the applicability of the results 
for the rates themselves. For example, this relation is independent of Xtc, the form and values of the 
supersymmetry breaking parameters and the form of the RGE above Mq (even if the theory becomes 
non-perturbative) . 

Contrary to the SO(IO) case, in SU(5), all the form factors in eq. (^) contribute in principle at the 
same general level. 

In practice, the contribution of the magnetic form factor g^" is almost always numerically dominant also 
in the SU(5) case, at least as long as the gaugino mass parameter M2 is not close to zero or the magnetic 
form factor has no accidental cancellation, which may occur for /i < 0. Relative to the penguin contribution 
gi, this comes about because the electric form factor has no term proportional to Af, and /xtan/3. The box 
contribution has a quite different structure from the penguin contribution. From the diagrams of figure |l^ 
one obtains for the effective Hamiltonian involving the quark fields 



ith 



cos' t/w 

1 4 1 

{cuL ■ ULl^UL + Cdx. • dL^dL) - 7:(c„R • urY^r) - -^{cdn ■ dB^^dn) 



(40) 



36 
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4 

E 

n.m—l 



[S(m?^, , M„, M„) - B(mi^, , Af„, M„)] x 

xH^BHmSiHnB + ^^^W, COt(?w)(i^™B + ^H,^w, COt ^w) 

[B (m?^ , mj^ , M„ , ) - B (m?^ , mj^ , M„ , M„ )] x 
xH^sH^BiHnB - 3H„w, cot0w)(i^™B - 3H,^w, cot^w) 

4 

J2 [B{ml^ , , M„, M„) - B{ml^ , , M„, M„,)]Hl^H^^ 



4 

E 

n.m—l 



n,m—l 
4 

E 

n,m— 1 



(41a) 

(41b) 
(41c) 

(41d) 



Consequently, by means of {N is the number of neutrons in the nucleus) 
one has 

'^^fi = a Jifl ^2 K;^;;^[^(32c„„ - 2c„, + 4cd, - CrfJ + iV(16c„, - + Sc^, - 2cdJ]. (44) 
(47r cos fw j I z 

A numerical calculation shows that the contribution of the box diagrams to the decay rate goes significantly 
below the analogous contribution from the penguins as soon as M2 moves away from zero, due to the rapid 
increase of the squark masses in the denominator and, even more so, to an increase of the magnetic 
contribution. In figure |l^ we give the rate for /i — > e conversion in the SU(5) case and tan (3 = 2. 
The numerical results for tan/3 =10 are not shown because they reproduce the relation ( |3^ ) to a good 
approximation for any value of the others parameters. 
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r 



/i7 



Very similar considerations to those developed in the previous sections can be made in the case of the 
T ^ decay. A main point is the relative difference between the SU(5) and the SO(IO) case. 

In the SU(5) case, the amplitude for r ^ fij is simply obtained from A^{ij. 67), equations ([1^19) 
with the replacement of the factor m^J/°*V^^ by mrV°*V°^, up to negligible terms of relative order m^/mr 
(remember that in SU(5) and /i^ are degenerate to a very good accuracy). Consequently the following 
relation holdaJ 



B.R.(t- 








B.R.(/x - 


67) 


SU(5) 


*^ er 



B.R.(t ^ivv) « 3 • 10^ 



0.77 

y 



0.01 



Vtd 



(45) 



with y given in (|5^). For given values of the mixing angles, this relation establishes the relative merit of 
the searches for the two decay processes as a possible signal of lepton flavour violation. With \Vtd\ = 0.01, 
the present limit on /i ^ 67 (B.R. < 4.9 • 10""'^^) |l^ is about 30 times better than the present bound on 
r — > /i7 (B.R. < 4.2 • 10^^) Q. Using the relation (p5|), the contours of figures ^ and || can be relabelled 
with values of B.R.(t /.i-f). 

In SO(IO), the ^ — > 67 amplitude is proportional to rrir and is therefore enhanced, as discussed in 
section 9. Consequently the ratio of the branching ratios will be further suppressed in SO(IO), relative to 
eq. (^5|), by an approximate factor of order (to^/tot-)^. 



12 Electric dipole moment of the electron 

It has been pointed out by Dimopoulos and one of us (L.H.) § that, in an SO(IO) unified theory, the low 
energy Lagrangian gives rise to an electric dipole moment for the neutron and the electron originating from 
the phases of the Yukawa couplings. We concentrate here on the dipole moment, de, of the electron, since 
in this case a very simple relation exists between de and the fj, ej rate. It is clear however that the 
search for a dipole moment of the neutron constitutes an independent and equally important signature for 
the general effect discussed in this paper. 

The full set of diagrams that contribute to the electric dipole or magnetic moments of the electron 
coincides with the one shown in fig. with A*l(a*l) replaced by eL(e^). In particular, as readily seen from 
the different dependence on the CKM matrix elements, only the diagrams of figures |l^b,c,c' contribute to 
the electric dipole moment (with V^^ replaced by V°e), since they are the only ones with an imaginary 
part. These are, on the other hand, the same diagrams that dominate the /i — *■ 67 amplitude through their 

^This equation corrects eq. (21) of ref. 



Figure 15: Same as in fig. ^except for the scale of the initial condition on the RGBs taken at 2.0-10^^ GeV. 



GUT 


r — !■ fij fi ^ e 




SU(5) 


0.03 ;^o.2 





SO(IO) 


< 0.03 0.2 


2siiitp 



Table 1: Relative merits of various observables relative to /i ^ ej. All branching ratios, as well as the value 
of de, are normalized to the present limits (B.R.(/^ 67) < 4.9-10"" [|l6|, C.R.(^ ^ e in Ti) < 10"^^ 0, 
B.R.(t ^ ^7) < 4.2 • 10-s |l|, 141 < 4.3 • 10-27 e . cm |l|). 



rrir dependent contribution. As a consequence, the following approximate relation holds between the form 
factor F2 defined in eq. (pi]) and the electron dipole moment, 



141 =e|F2| 
with the CP violating phase (p defined by 



sinif — e\F2 



td 



sm If 



(46) 



ImK(V;°e)'(KTn ^ |m.(K'=e)'(K°;)'|sin^. 
In SO(IO) the electric dipole moment is therefore approximately related to the /i — > 67 branching ratio by 

141 



10-27, 



1.3 sin ip 



B.R.(^ 67) 
1F12 ■ 



(47) 



Using this relation, the contours of figures |ll| and |l| can be relabelled with values of |4|/siniy9. It 
is interesting to notice that the present upper bound on /i ^ 67 (B.R. < 4.9 • 10"^^) and 4 (|4| < 
4.3 • 10-2^ e • cm) pi are almost exactly equivalent for shup = 1/2. 

As in SO(IO), in SU(5) too, the diagrams that could contribute to the electric dipole moment of the 
electron are obtained from those of figure ^ by replacing the muon with the electron in the external line. 
This time, however, no electric dipole moment arises since the CP violating phase disappears from the 
product V°*V^^ of the relevant CKM matrix elements. 



13 Conclusions 

In this paper we have discussed the lepton fiavour violating processes and electric dipole moments induced 
in a supersymmetric unified theory by the large top Yukawa coupling. Under the stated assumptions, the 
experimental study of these processes provides a very significant test of supersymmetric unification. Already 
the present experimental limits give, especially in the SO(IO) case, significant restrictions on the allowed 
parameter space, often considerably stronger than those inferred from direct searches of supersymmetric 
particles. 

The main results of this paper are the contour plots for B.R.(/i — > 67) shown in figures |^ and |^, for 
SU(5), and in figures 0and^, for SO(IO). The figures0and| can also be used for t ^7 by a relabelling 
of the contours using equation ( ^ ) . Similarly the contours of figures |ll]and|l| can be relabelled using ( pq ) 
and (Erh so that they apply to /i ^ e conversion and to de respectively. The case oi fi ^ e conversion in 



SU(5) is shown in figure 13. These plots cover the entire physical ranges of the parameters Ae and M2, 
and show the behaviour for both signs of /i and for both large and small values of tan/3. For large mg^ 
the B.R. decreases as l/m| ; however values of mg^^ above 400 GeV require a significant amount of fine 
tuning |2^. 

The sensitivities of the various processes to the SU(5) and SO(IO) theories is summarized in table |^, 
relative to that of /i ^ 67, with all observables normalized to the present experimental bounds. For 
sint^ < 1/2 all entries of this table are less than unity, showing that, for this case, /j, — > 67 is presently the 
most powerful probe in all cases. For simp > 1/2 the electron electric dipole moment provides the best 
probe of the SO(IO) theory. The decay t fij will only become competitive with the construction of a 



T factory. Future technologies and experimental possibilities should allow an interesting competition to 
develop amongst the other three processes. 

An additional prediction of this work is that the mass of f/j is suppressed significantly beneath that of 
jlii and Cfl, as can be seen in figure ^. In addition, in SO(IO) the mass of is suppressed beneath that 
of /2l and e^. This result is important for superpartner searches at e'^ colliders: the lightest charged 
scalar superpartner is almost certainly a scalar tau. 

The sizes of these effects depend on the following two main assumptions: 

i) The value of the top Yukawa coupling at the unification scale is large; 

ii) The field theoretic renormalization group equations, valid at the unification scale, can be extrapolated 
without substantial modifications up to the reduced Planck mass, Mpi = 2.4 • 10^"^ GeV, or to a scale 
close to it, e.g. the compactification scale of string theory. 

A value of the top Yukawa coupling at the unification scale less than one leads to a substantial reduction of 
the rates, as shown, e.g., in fig. ^as compared to fig. Also significant, although relatively less important, 
is the lowering of the scale for the universal conditions on the supersymmetry breaking parameters, as 
exemplified in figure |l5| where such scale is taken at 2.0 • 10^^ GeV. 

What other features of the unified model influence our results? Other than on the gauge group itself, 
the lepton flavour violation effects certainly depend on the specific form of the flavour interactions. In 
this paper we have studied the two simplest unified sectors that we know, one in SU(5) and the other 
in SO (10). Each has the minimal number of flavour Yukawa matrices, A" giving mass to up quarks and 
A'* to down quarks and charged leptons. It is well known that these Yukawa couplings lead to unrealistic 
mass relations between the light fermions. However, as already discussed in ref. ||], we do not expect that 
the necessary modifications of these couplings may lead to significant suppressions of the lepton flavour 
violation processes. They will rather give rise to an increased range of predictions about the central values 
discussed here. 

The traditional probes of supersymmetric unified theories are provided by proton decay, neutrino masses 
and by predictions for quark and charged lepton masses and mixings. These probes also have model 
dependences which arise from the choice of the gauge group and the flavour interactions, as discussed 
above. However, for each of these three probes, there is also a much greater uncertainty than in the lepton 
flavour violation processes. A generic unifled model has a free parameter for each of the flavour masses 
and mixing parameters of the standard model, and hence does not make predictions for the quark and 
charged lepton masses. Such predictions only arise when the form of the flavour interactions are restricted 
by further assumptions. 

While proton decay and neutrino masses are generally to be expected in superunified models, the sizes 
of these signatures are extremely model dependent. Consider first the case of proton decay. All superunified 
models contain baryon and lepton number violating interactions which couple the quarks and leptons to 
a set of superheavy coloured states H. The amplitude for proton decay depends on the mass matrix for 
these H states. This is perhaps the least understood, and most model dependent, feature of superunified 
theories, because it is directly related to the problem of why the Higgs doublets are much lighter than 
Mq. Only in one particular model |^, where the Higgs are made light by an extreme fine tune, has it 
been possible to relate the H mass to known parameters of the theory and hence make predictions for the 
proton decay rate. In fact the resulting rate is large, and this model is close to being excluded. In many 
other models the matrix structure of the masses for the H states leads to a large suppression of the proton 
decay amplitude, which then becomes gauge dominated, yielding a rate which is expected to be about four 
orders of magnitude below present experimental limits. 

The three neutrinos frequently acquire small masses in superunified models, particularly if the gauge 
group contains SO(IO). However, the size of these masses is inversely proportional to M^, the Majorana 
mass matrix for the right-handed neutrinos, which breaks lepton number and is typically not directly related 
to known parameters of the theory. A simple expectation of Mn « Mq gives masses for i/g, v^, which 
are too small to see in accelerator or reactor experiments. 

By comparison with these great uncertainties, which afflict the traditional signatures for superunifled 
models, the model dependence of the rates for Lg, L^, Lr and CP violating processes discussed here seems 
quite mild. 

We therefore conclude that searches for the Li and CP violating signatures discussed in this paper 
provide the most powerful known probes of supersymmetric quark-lepton unification with supersymmetry 



6^ = {3, 3, 9} c" = 96/5 c« = 24/5 

^u^oj fed ^{0,0,3} €^^ = 84/5 c^ = 36/5 

qnnm ^-o- {4, 4,14} c" = 63/2 c* = 9 

s'JUu; fed = {0,0, 10} c<i = 63/2 c*=45/4 



Table 2: Values of the RGE coefficients in SU(5) and in SO(IO). 



breaking generated at the Planck scale. For example, an experiment with a sensitivity of 10~^^ to B.R.(/i — > 
e"f) would probe (apart from a small region of parameter space where cancellations in the amplitude occur) 
the SU(5) model to A(g = 1-4 and mg^ = 100 GeV, and would explore a significant portion of parameters 
space for mg^ = 300 GeV. In the SO(IO) case, where the present bound on — > 67 is already more stringent 
than the limits from high energy accelerator experiments, a sensitivity of 10~^^ would probe the theory to 
A(G = 1-25 and nig^ close to 1 TeV. 
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A Renormalization from the Planck to the GUT scale 

Neglecting all couplings except the gauge and the top Yukawa ones, the solutions to all the one loop RGBs 
between i^max — My>\ and i?min = Mq can be given analytically. 

The RGBs for the dimensionless couplings and for the dimension-one soft terms are 

(48a) 
(48b) 
(48c) 
(48d) 

where t{E) = (47r)^^ lnMpj/£'^, 55 is the coupling constant of the unification group, Q!5 = 5|/47r, M5 is 
the gaugino mass, g = 1, 2, 3 is the generation number and the values of the numerical coefficients in SU(5) 
and in SO(IO) are given in table ^. The subscript '5' stands for 'unified' rather than for SU(5). We also 
set &3 = 64. 

The full analytic solutions of these equations with the boundary conditions 

a5(MG) = ac, Mz{Mg) = Mgc, ^(Mg) = A^g 

and universal A-terms at the Planck scale, Aq, are 

a^{E) = f^\E)-aG (49a) 
AME) = f^\E)-M,G (49b) 

\ 2max f Tp\ 

>h{E) = ^ ^-^ -r (49c) 

1 + A^-(£;)(ArG' - X^t'mr ' ""{E) 

AliE) = Ao + xl{E)M5G-hp'{E)/ht (49d) 
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A'iiE) = Ao + xi{E)M5G-b''j'{E)/bt (49e) 



where the functions fz{E), x^{E), X"^'^^{E), I'{E) are explicitly defined below in eq. (50). 

Assuming universal values at the Planck scale for the dimension- two supersymmetry breaking soft terms, 
m|j = m§, the one loop RGBs in SU(5) 

|™|^ = 2c"9lMl 

= 2c^glMi-^jxU2m%+mj, + A',) 



dt 



are solved by 



m\{E) = 
m]j{E) = 



m 



mUE) = ml + xliE)Mic-IiE)Sg3 



and in SO(IO) 



dt 
d 

— IT, 

dt 



4 
14 
5 

14 



by 



4{E)Ml^^2-I{E) 



,^,{E)Ml^-i-I{E)5,, 



In both cases, we have defined 

h{E) = l + glbG[t[E)-t{MG)] 



gliE) 



1 - [/5(i?max)//5(i?)]l+^"/^'= 



J(i?) = p 



i(l-pK-^(l-p)(l-M: 



2max 
t 



t)AoAk 



5P1 



I'{E) 



p[Ao - (1 - btXf''^H)M, 



5P1J 



(50a) 
(50b) 

(50c) 

(50d) 

(50e) 



and p{E) = \^{E)/ Xf^'^^^E) < 1. From eq.s (50), one learns that the main factor that determines the 
size of the lepton flavour breaking parameters / and /' is the overall factor p. In turn, p is only weakly 
dependent on the /?- function coefficient (see figure ^) . 

For the numerical values at AIq of the different quantities defined above we take 

ac = 1/24, Mg = 2.0 • 10^^ GeV, Afpi = 2.4 • 10^^ GeV, 6g = -3, 

so that tQ — 0.0606 (a subscript 'G' on the various functions of E indicates that they are evaluated at 
Mg). In SU(5) AJ^^'^ = 1.56 and 



Ig = /(Mg) 
EE I'{Mo) 
while in SO(IO) Xf}?' 

Ig EE /(Mg) 
I'g EE I'{Mg) 



PG 



\{l - pg)AI + 0.198(1 - pg)^oA/5G + (0.224 - 0.029pg)M|g1 (51a) 



= pG[Ao + 0.298Af5G] 
~ 1.36 and 



PG 



ml + -(1 - pg)AI + 0.343(1 - pg)^oM5G + (0.435 - 0.088pg)M5G 



/9g[^ + 0.515M5g] 



(51b) 

(52a) 
(52b) 



B Renormalization in the MSSM 

Neglecting all couplings except the gauge and the top Yukawa ones, the solutions to all the one loop RGBs 
between -Emax = Mg and i^min = Mz may be written in terms of analytic functions and only one function, 
Xf'^^{E), calculable only numerically pl| . 
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Table 3: Values of the RGE coefficients in the MSSM. 



The RGEs for the dimensionless couplings are 

d 1 , , 

-— - 4Trh 53a 
dt Ui 

= A?(q"5f "6tA?) (53b) 

where i = 1,2,3 runs over the three factor in the Standard Model gauge group U(l) (E) SU(2) (g) SU(3), 
t{E) = (47r)^^ IuMq/E'^ and the values of the coefficients are shown in table ^. 
The solutions with boundary conditions ^^(A/q) = Q!q and Xi^Mq) = Ate £^re 



a,iE) = fr\E)-aG (54a) 

;)^2max(-^) 



\ 2max f T7i\ 

>h{E) = }J (54b) 



where ii{E) = 1 + h,g%t(E) and 

/•InMc 

£;„(£;) = 17/^' (^): ^u(i?) = 2 / E,,{E)dhiE, 

J\n E 



The Yukawa couplings of the fcrmions in the diagonal basis scale as Xg{E) — Xg{MQ) ■ y^'s eI/^ where 
g = 1, 2, 3 is the generation number, a = u, d, e and 



y{E) = exp 



In Mo \2( pi 



\nE 167r^ 



[l-p(i?)r-s p(i^) ^ < 1. (56) 



The factor y used in the text is given hy y = y{Mz)- 

The RGEs for the three gaugino masses M^, the supersymmetric /i-term and the A terms are 



d Mi 
dt ai 


= 


(57a) 


d 

dt^ 


= i(2c,^,2_fou^2)^ 


(57b) 


^A 

dt '''' 




(57c) 


d 

Jt^"'' 


= 4gfM,-b^gX^,Ad,g 


(57d) 


'^A 


= Cig^M, 


(57e) 



nth all the various coefficients listed in table The solutions are 

M,iE) ^ fr\E)-M^G (58a) 

^i[E) - M(M5G)-2/'"(S)Sh(^) (58b) 

A^iE) = Al^+xUE)M,G~by{E)/bt (58c) 

Al[E) = A%+x'i{E)M^G-bp\E)/bt (58d) 

Al{E) = A^G+x?(ii;)Af5G (58e) 



with x^{E) and /'(£') defined below in (60) and Eh{E) in (||). 

The RGEs for the dimension- two soft parameters of a representation R — {Q, u, d,e,L — h} are 



dt ^ 



except for the multiplcts ft-u, Qs and t involved in the top Yukawa coupling. For them 



dt 
d_ 

d 2 
—mr 
dt * 



= 2ctgfMf-^bt\UA't+3m') 



where m'^{E) = [m\^{E) + to^ (S) + m~(i?)]/3. The solutions are 

m2,(i?) - m|,,(AfG) + X2«M|G 
except when the top Yukawa coupling appears, where 



m%(E) 



mUE) 



mliMo) + x^^{E)Ml^ - \iI{E) 



'103 
2 



(Mg) 



5G 
2 



3/(i?) 



mi^^{Mo)+xl{E)Mi^~-mE) 



(59a) 

(59b) 
(59c) 
(59d) 

(59e) 

(59 f) 
(59g) 
(59h) 



and 



x^{E) 

m 

I'{E) 



E7?[/r"(^ma.)-/r"(i^)] 



m^{Mo) + \{l - p)Alo \{l p){l h\ 



2max 
t 



t)M5GAtG- 



1, 



[p(l - btAr^'^i)' - h\r^H\c^gi)]M, 



5G 



p[^tG - (1 - btXt^-^t)M, 



5GJ 



(60a) 
(60b) 

(60c) 



Notice that, apart form obvious replacements, a;^, / and /' maintain exactly the same form as in eq. (50). 

The numerical values at Mz of the different quantities defined above are t{Mz) = tz = 0.418 (a 
subscript 'Z' on the various functions of E indicates that they are evaluated at Mz), 



E„ 



13.6, btF^z = 10.5, 



\ max 



1.14 



Iz = pzm{M^) + ^{l-pz)AjQ + 1.50pzil-pz)M5GAtG + Pzi4..37-1.70pz)Mlc (61a) 

I'z = pz[AtG + 2.26M5G] (61b) 

/i2 = 2.44 /2z = 1.22 /32 = 0.343 
x'lz = 4.02 xfz = 3.98 x^^ = 0.700 
x^z = 7.16 x^z = 6-73 x^z = 6.68 x^z = 0.151 x^z = 4z = 0.528 

Finally the Higgs doublets mass parameters = m\^{Mz) and = (M^) defined in eq. (10) may 
be expressed in terms of the universal supersymmetry breaking parameters as 

ml=ml + (xfc + x'izWlc, ml = ml ~ ^Iz + Ig), (62a) 

in SU(5), while, in SO(IO) 

ml^ml + (x|g + x'^z)MiG, ^1 = ml - {^z + f/c). (62b) 
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